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Abstract 

We consider a two-component Bose-Einstein condensate in a one-dimensional optical cavity. 
Specifically, the condensate atoms are taken to be in two degenerate modes due to their internal 
hyperfine spin degrees of freedom and they are coupled to the cavity field and an external 
transverse laser field in a Raman scheme. A parallel laser is also exciting the cavity mode. When 
the pump laser is far detuned from its resonance atomic transition frequency, an effective nonlinear 
optical model of the cavity-condensate system is developed under Discrete Mode Approximation 
(DMA), while matter-field coupling has been considered beyond the Rotating Wave Approxima- 
tion. By analytical and numerical solutions of the nonlinear dynamical equations, we examine 
the mean cavity field and population difference (magnetization) of the condensate modes. The 
stationary solutions of both the mean cavity field and normalized magnetization demonstrate 
bistable behavior under certain conditions for the laser pump intensity and matter-field coupling 
strength. 
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I. INTRODUCTION 



Recently, due to experimental advances in coupling a dilute gas of bosons to a single 
mode of an optical cavity, many theoretical and experimental works have been performed 
in order to explore and explain the physics of such complex systems. Bose-Einstein con- 



densate (BEC) itself is a very rich platform which a 
of quantum systems such as quantum turbulence 



lows us to examine different properties 

in n 

1H3(, quantum chaos [4j and entangle- 



ment [5|. A Bose-Einstein condensate can also be used to study nonequilibrium dynamics 
and decoherence in finite quantum systems [6j as well as calculation of fluctuation indices for 

n i 

atomic systems [7|. Moreover, if confined in an optical lattice [8], Bose-Einstein condensates 



10], 
HI, 



provide the possibility to observe solid-state physics processes, such as localization 
in these atomic systems. Further uses of BECs are exemplified in transport problems 
interaction driven instabilities {12I. stability of boson- fermion gaseous mixtures [13I] , and 
Raman pump-probe experiments [14]. 

In a more complex setup, where a dilute condensate of bosons is confined inside a high- 
finesse optical cavity subject to external laser fields, if the cavity mode and laser fields 
are far-detuned from the atomic transition frequency of the condensate atoms, system is in 
dispersive regime. Under the dispersive regime conditions, atom-field interaction provides an 
optical lattice for the condensate atoms which affects their mechanical motion. On the other 
hand, the atoms cause a position-dependent phase-shift of the cavity mode. As a result, 
the condensate-cavity system is highly nonlinear with nonlocal nonlinearities which give rise 
to a series of interesting physical phenomena such as self-organization of condensate atoms, 



m 



161 ] and optical 



cavity enhanced superradiant scattering, Dicke quantum phase transition 
bistability. The nonlinearity caused by atom-filed interaction can be more effective (17| than 
the nonlinearity from atom-atom interaction, which itself can play crucial role in scattering 
dynamics of condensates . 

The optical bistability, which is the focus of this article, has been studied in spinor BECs 
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201 ] and two-component BEC with two modes coupled by a classical field (2l||. In a 
system consisting of single-mode BEC in an optical cavity, a transverse laser pump has also 
been used to control the bistability of cavity photons induced by a parallel pump [l7j . In this 
work we consider a two-mode BEC in a one-dimensional optical cavity, where two laser fields 
are applied to the system, one parallel to the cavity axis and the other one perpendicular 



to it. Specifically we assume that the transverse pump is scattered to the cavity mode by 



condensate atoms in the Raman scheme 
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23] . Under this condition, we examine the 



effect of transverse pump strength on the bistability of both mean cavity photon number 
and normalized population difference (magnetization) of the two modes. Stationary state of 
cavity field and condensate wavefunction are obtained under Discrete Mode Approximation 
(DMA) which has been shown to be reliable for similar systems 17|, \24 \. 

This paper is organized as follows. In SecJTTI first the model for our system is introduced, 
then the Hamiltonian of the system and equations of motion of cavity field and condensate 
are derived. Using DMA, we solve the equations of motion for the steady state of the system 
in Sec JIIIt where we show how mutual bistability of the mean cavity photon number and 
magnetization take place under certain conditions for laser field intensities and cavity-atom 
coupling strength. Finally, we summarize our work in Sec JIVI 



II. TWO-MODE BEC IN ONE DIMENSIONAL CAVITY 



We consider a condensate of N atoms, each with two internal degrees of freedom, shown 
as states b and c and an excited state e, in a one dimensional cavity along the x axis. The 
cavity has a single mode with frequency u c and is subjected to laser field with frequency 
Uq, which is far detuned from the atomic transition, in directions parallel and perpendicular 
to its axis. If the transverse laser field interacts with matter through Raman scattering 
(Fig. 1), then the full Hamiltonian of the system will have following form 



H = 2 j dx $ + v ^ + hwbc6 ^) ^ + Yl J dx 

j=b,c i,j=b,c 

+ hw^a - ih m {ae iulot - aV^ ') + H Raman , (1) 



where iOf, c is the frequency of transition between modes b and c and are the interaction 
strengths of atoms in modes i and j. The parallel laser field intensity is shown by 7711 and a 
and a) are the annihilation and creation operators of the cavity mode. The interaction of 
atoms with the transverse pump is shown by the Raman scattering Hamiltonian (HRaman) 
and has the following form 
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FIG. 1. Left: Schematic drawing of a BEC in one dimensional optical cavity subject to parallel 
and transverse laser fields. Cavity has a decay rate of k. Right: Two internal modes (b and c) 
of the BEC atoms are coupled by cavity field and laser field via the atomic excited state e in 
Raman scattering manner. Both laser field and cavity field are detuned from the atomic transition 
frequency. 



H Ramm = -to J dx 4,%(e-"> 1 + e'""'}4, b + H.c. 

- ifi / d ^I 90 cos( fa )(« + ^ + tf,., (2) 

where ho and go are the atom-pump and atom-cavity coupling strengths, respectively, and 
the Hamiltonian is written without Rotating Wave Approximation. After adiabatically 
eliminating the excited state ip e and introducing U = go/A and r\ = h g /A , with A = 
u — Ube being the pump detuning from the atomic transition, one obtains 



H Raman = h V jdx(a + at)(e" iwo * + e^ ot ) cos(kx) + ^ 

+ f dx ( e — ot + 



A 



+ hU J dx cos 2 (kx)(a + a j ) 2 ^pltp c . (3) 

Substituting ([3]) into (PQ), moving to a rotating frame defined by the unitary operator 
U = e~* tJ ° iata , and ignoring two-photon processes result in 

H = E J dx + E / dx ^M-^M-. ( 4 ) 



j=b,c i,j=b,c 
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for the Hamiltonian of the system where 



h 2 d 2 

H = 77— r + ( frUo cos 2 (kx)(aa^ + a) a) + VJx) + hujbc) <y + cr 

2m ox 2 v ' 

flhh 2 \ 
+ — — - + V b (x) a~a + + fvq{a + a ] ) cos(kx)(a~ + a + ) 



. A 

— HS c a/a — iHrf\\ (a — a'). (5) 

Here a + = \ip c x ^| and a~ = \tpb x ip c \ are the ascending and descending operators in the 
two-mode manifold and 5 C = uo — u c is the pump-cavity detuning. From equations (jlj) and 
(jSJ) one can derive the Heisenberg equations of motion for condensate modes and cavity field 
as follow 



; i ( h 2 d 2 . 2hh 2 , t , (t , \ 

A = - ^ ( -2^q-2 + V *( X ) + ^ + U »W + U&<»c 1 

— irj cos(kx) (a + a^)tjj c , (6) 



if h 2 d 2 hg 2 , \ 

ip c = -- ( - — + V^z) + ftu; fec + cos 2 (/cx)(aa t + a f a) J ^ c 

- ^ \ Ucciplipc + Ubctplipbj *Pc - if] cos(kx) (a + a 1 )^, (7) 

and 

a = i (in + 5 C ^ / <ir ^ cos 2 (kx)ib c I a 

where we have introduced a decay rate k for the cavity in (jSJ). In order to obtain the 
equations of motion for the expectation values of field operators, we further simplify our 
system. Since we work with a condensate of large number of atoms, we can safely treat 
the condensate as a coherent state and disentangle it from photons. We also neglect higher 
order correlations within the photon subsystem and use < a >< a> > and < >< a > 
instead of < aa^ > and < a^a >. At the end, as far as we only deal with the intensities, we 



treat the photon as a coherent state as well. As a result, following equations of motion are 
obtained 



% ( & 32 x r ( \ 2hh l , / 12 1 , ,2^ , 

^ b = -h\-2md^ + Vb{x) + ^ + Ubb ^ +UM 

— 2ir\ cos(kx)a r ip c , (9) 



% ( h 2 d 2 \ 
ip c = - - + v c( x ) + fi^bc + 2^f/ cos 2 (A;x)|a| 2 j ip c 

- t (u C c\4>c\ 2 + u bc \?p b \ 2 ) i) c - 2ir] cos(kx)a r ip b , (10) 



a = i yiK + 5 C — 2Uq J dx \ijj c \ 2 cos 2 (kx) J a 

-irjjdx co8(kx)(ip2i>b + iPt^c) + V\\, (H) 
where a r is the real part of the cavity field a. 



III. BISTABILITY OF PHOTON NUMBER AND MAGNETIZATION 

To solve the equations of motion derived in SecJTT]for the steady state of the system, we 
employ Discrete Mode Approximation. The ground state of the condensate, without any 
laser field, is a homogeneous macroscopic state with zero momentum, which we refer to by 
0o- Above superradiance threshold one can assume that the condensate is fragmented to a 
symmetric superposition of the states with ±hk momentum due to the transverse laser field. 
On the other hand, since the cavity mode is excited by the parallel pump, absorption and 
emission of the cavity photons can excite the condensate to the superposition of states with 
±2hk momentum. Therefore if we consider first order perturbation on the homogeneous 
wavefunction <p , following functions can be used as the basis for DMA: 



0o= Vl/L 



>i = a/2/L cos(kx) 
; = ^/2/Lcos{2kx) 



(12) 



with L being the length of the cavity (condensate). Now the wavefunctions of the two modes 
of the condensate can be expanded in this basis as follow 



ip b (x,t) = ^2<j>ibi 

i=0 
2 



(13) 



If we substitute these wavefunctions into the equations of motion of the two modes, 
by ignoring external potentials and atom-atom interaction, the equation of motion for the 
condensate wavefunction can be written in a compact form: 



ih^-X = H(a)X = [H + H 1 + H 2 + 2aH 3 + 2a r H i }X, 
at 



(14) 



where X = (6 , &i, &2, Co, ci, c<i) T and 



Hn = hu r 
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Here u r = hk 2 /2m is the recoil frequency. 

To examine the equilibrium properties of the system, we set a = in ffTT]) which results 

in 



a 



T] J dx COs(/cx)(^*Vb + TpblPc) + ^|| 



(20) 



in + 5 C — 2U J dx cos 2 (kx)ipcipc 
On the other hand, one can easily obtain the following alternative expressions for the inte 
grals in the above equation 



dx cos 2 \kx)il)*ip c 



hu ' 



(21) 



S 



/ 



dx cos(kx)(ip b *ip c + ip* c ip b ) = — X f F 4 X 



(22) 



Therefore cavity field a and averaged photon number n = \a 



2 will read as 



\X^H A X + ir)\ 



(23) 



a = 



n = 



^H.Xf + r, 2 



(24) 



^ + {s c -^xm 3 xy 



Here one should notice that if there is no transverse laser field (ho = 0) then rj = and 
H4 = 0. As a result, there will be no transition between the two modes and if atoms are 
initially in the ground state b, they will stay there. Therefore expectation values of both 
i?4 and H 3 in ( 124"]) or, equivalently, the integrals in ( 120]) would be zero and no bistability 
in the system is expected. The importance of the transverse laser pump could be predicted 
from the crucial role it plays in Raman scattering Hamiltonian. On the other hand to 
better understand the contribution of parallel laser field in bistability, we remind that due 
to dependence of X on a, expectation values in numerator and denominator of (l24"j) are 
functions of n. To lowest order, we can assume that these terms depend on n linearly. 
Therefore ( 124"]) is a cubic equation of n which in the case of r]\\ = has always a zero root. 
Therefore, to avoid a zero root in (]24]) . we will consider cases where r]\\ ^ 0. 

To obtain the mean cavity field n for the steady state of the system we need to find the 
wavefunction X s in steady state. Since H(a)X s = E$X S is nonlinear due to the dependence 
of X on a, we first solve it for Eq and X s with a guess for the value of a. Then by substituting 
the resulting X s into (]23]) . a new a is obtained. If this new a is equal to the guessed value 
then the steady state is reached, otherwise we repeat the procedure, while using this new a 
as the guess, until the steady state is attained. 

In addition to the average photon number, in a system consisting of two-mode condensate, 
there is another quantity which reveals the nonlinear effects of matter-field interaction and 
that is the normalized population difference (magnetization Z) of two modes. In our system 
magnetization is defined as Z = J dx(\ipb\ 2 — \ip c \ 2 )/N, where N = J dx(\ipb\ 2 + \i^ c \ 2 ) is the 
total number of atoms and is fixed. 
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FIG. 2. Mean cavity field n (left) and normalized magnetization Z (right) as functions of 77 for 
different values of parallel pump strength 7711. We remind that 77 is proportional to transverse 
pump strength and atom-cavity coupling strength. In all plots Uq = —0.5 u r , Aq = —4 x 10 6 oj r , 
k = 400 u r , 8 C = 4800 oj r , N = 4.8 x 10 4 , and oj bc = uj t . 

Fig. 2 shows mean cavity field n and normalized magnetization Z as functions of 77 for 
different values of parallel pump strength rj\\. We have considered a condensate of N = 48000 
atoms in a cavity with decay rate k = 400 u r . The laser field is detuned from cavity mode 
by S c = 4800 u r and from atomic transition by A = —4 x 10 6 u r . In all plots, atom-cavity 
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— T|/C0 ~ r |/<» r 

FIG. 3. Mean cavity field n (left) and normalized magnetization Z (right) as functions of 77 for 
different values of atom-cavity coupling strength Uq. In all plots 7711 = 500 oj r , Aq = —4 x 10 6 u r , 
k = 400 oj r , 5 C = 4800 u r , N = 4.8 x 10 4 , and oj bc = un- 
coupling is assumed to be U = —0.5 co r . As one can observe from these plots, the effect of 
parallel pump on widening the area of bistability is not monotonic. While increasing 7711 from 
10 to 500 results in a wider interval of rj in which bistability happens, further increasing it 
to 1000 has the opposite effect. Moreover, 7711 has larger impact on the width of the area in 
which bistability takes place than on the values of n and Z in bistable region. The interesting 
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point about this system is the fact that nonlinear effects of matter-field coupling result in 
bistable behavior of both matter and field, such that, bistability happens for magnetization 
of condensate atoms exactly at the same region of rj where n has shown bistable behavior. 
Moreover, it can be seen in Fig. 2 that, in some points, two states not only with different 
values of magnetization but also with different signs of magnetization are stable. 

In the next step, to better understand the role of atom-cavity coupling on bistable behav- 
ior of the system, we keep the value of parallel pump strength r)\\ constant and increase the 
strength of atom-cavity coupling. Fig. 3 shows the average photon number n and normalized 
magnetization Z as functions of rj for the same values of parameters used in Fig. 2 but this 
time with 7]\\ = 500 ou r and for three different values of Uq. Clearly, atom-cavity coupling has 
a significant effect on the width of the area in which bistability happens. More importantly, 
Uq can change the distance between the two stable branches. It is worth mentioning that r\ 
itself is proportional to go and therefore for a fixed value of 77, larger Uq means smaller value 
for h . As a result, increasing U makes it possible to achieve bistability with smaller values 
and wider range of atom-laser coupling h . 

IV. CONCLUSION 

In this work we developed an effective Hamiltonian and equations of motion for a cavity- 
condensate system consisting of a two-mode BEC in a one dimensional cavity, while parallel 
and transverse laser fields are applied to the system. Under DMA, simultaneous and mutual 
bistability (multistability) of cavity field and population difference (magnetization) of the 
two modes has been observed for different values of transverse and parallel pump strength. 
The system shows bistable behavior for quite wide range of parameters. Moreover it has 
been shown that, with strong enough cavity-matter coupling strength, bistability happens 
between two states with different sign of magnetization. 
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